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Abstract 

Many geophysical processes can be modelled by using interconnected networks. The small- 
world network model has recently attracted much attention in physics and applied sciences. In 
this paper, we try to use and modify the small-world theory to model geophysical processes such 
as diffusion and transport in disordered porous rocks. We develop an analytical approach as well 
as numerical simulations to try to characterize the pollutant transport and percolation proper- 
ties of small-world networks. The analytical expression of system saturation time and fractal 
dimension of small- world networks are given and thus compared with numerical simulations. 
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Introduction 

Many geophysical processes such as flow in disordered porous media, magma transport and 
saltwater intrusion [Holzhercher,l'd'd^] can be modelled as percolation and diffusion in inter- 
connected networks that can generally be represented as connected graphs. The properties of 
complicated networks are mainly determined by the way of connections between occupied sites 
or vertices. Two commonly used quantities, which characterise a graph, are the degree of local 
clustering and the average distance between vertices. The generic structure and the degree of 
the randomness of networks have strong influence on the overal properties. Obviously, regular 
network is the one limiting case with a high degree of local clustering and a very large average 
distance, while the other limiting CEise is the random network which shows negligible local clus- 
tering and the average distance is relative small. In reality, most networks are neither exactly 
regular nor completely random, realistic networks are something between these two extremes. 

Small-world network is a special kind of networks with a higher degree of clustering and a 
small average distance. The degree of the randomness of networks often has a strong influence 
on the network properties. Recently, Watts and Strogatz [1998] has presented a mathematical 
model to characterises the network structure and to study the dynamic properties of small- 
world networks. Such small-world phenomenon can be obtained by adding randomly only a 
small fraction of the connections, and some common networks such as power grids, fllm stars 
networks and neural networks behaves like small-world networks [Newman and Watts, 1999; 
Moukarzel,1999; Newmann et al, 2000, Boots and Sasaki, 1999; Barthelemy and Amaral, 1999]. 
Most recent studies focus on the spreading and shortest paths in system with sparse long-range 
connections by using the small- world model [Newman and Watts, 1999; Moukarzel,1999; Pandit 
and Amritkar, 1999] . The spreading of some influence such as a forest fire, an infectious disease 
or a particle in percolating media is studied using a simple rule: at each time step, the influence 
propagates from the infected site to all uninfected sites connected to it via a link although this 
link is not necessary physical. A long-range connection or shortcut can simulate the spark that 
starts a new flre spot, the infect site (say person with flu) suddenly travel to a new place or 
site, or a portable computer with virus that start to connect to the network a new place. These 
research mainly focus on the immediate response of the network by adding randomly some long- 
distance shortcuts, there is no time delay in the network systems. However, in reality, a spark 
or an infection can not start a new fire spot or new infection immediately, it usually takes some 
time A, called ignition time or waiting time, to start a new fire or infection. Thus the existing 
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models arc no longer be able to predict the response in the networks or systems with time delay. 
Yang [2000] has studied the effect of time delay on the response of small-world networks. 

The purpose of the present work is to model geophysical processes such as pollutant trans- 
port, permeability model and saltwater intrusion [Holzbercher, 1998] as a particle percolation 
process on the small-world networks with time-delay by using or modifying the small- world the- 
ory so as to characterise the small-world networks with time-delay by using fractal dimension 
and system saturation time in the percolative media. The present model will generally lead to 
a delay differential equation, whose solution is usually very difficult to obtain if it is not impos- 
sible. Thus the numerical simulation becomes essential. However, we will take the analytical 
analysis as far as possible and compare with the results from numerical simulations. 



Mathematical Model of Small- World Networks 

Considering a randomly connected network with n-dimensional lattice [Watts and Strogatz, 1998; 
Bollobas, 1985; Watts, 1999] and n = 1,2,3. Overlapping on the ra-dimensional networks are 
a number of long-range shortcuts randomly connecting some vertices, and the fraction of the 
long-range shortcuts or probability p is relative small p <C 1. Now assuming an influence or a 
pollutant particle spreads with a constant velocity w = 1 in all directions and a newly infected 
site in the other end of a shortcut will start but with a time delay A. Following the method 
developed by Newman and Watts [1999] and Moukarzel [1999], the total influenced volume V{t) 
comes from two contributions: one is the influenced volume with Vd t^~^dt where t is time, 
and the other contribution is 2pV{t — C, — 1\) for a sphere started at time C,. By using a continuum 
approach to the network, then V{t) satisfies the following equation with time delay 

v{t) = vj r-^[i+r'V(f-c-A)]dc, (1) 

Jo 

where n = 1, 2, 3 and is shape factor of a hypersphere in n-dimensions. The length scale of 
Newman- Watts type [1999] can be defined as 

where k = const being some fixed range. By proper rescaling t by (Fdw(n — l)!)^/", that is 

r = t(Fdr"(n-l)!)'''", 5 = A(Fdr"(« - I)!)'''"- (3) 

By rewriting (1) in the rescaled form and differentiating the equation n times, we have a time- 
delay equation 

'£Y.=C + V{r-S), (4) 

which is a delay differential equation, whose explicit solutions is not always possible depending 
on the initial conditions. For a proper initial condition V{t) = for — 5 < t < 0, the solution 
can be explicitly written by using the time-stepping method, 

t/<5 

Clearly, for t ^ and 5—^0, the our solution degenerates into 

oc 

given by Moukarzel [3]. The solutions explicitly become V{t) = ^(e^ — 1) for n = 1, and 
V{t) — 5^ (cosh r — 1) for n = 2. However, when 5 7^ 0, there is no such simple expressions. 
It can be expected that the time-delay parameter 6 will have a strong effect on the evolution 
behaviour V{t) of an influence. 
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Pollutant Transport and Saturation Time 



By using the small-world network model, wc can sinmlatc the diffusion and transport process 
in the disordered media. Imagine a pollutant particle walk randomly on the networks which 
correspond to the interconnected networks in a way structured like the swiss-cheese model for 
disordered porous media, while the longe range shortcuts corresponds to the fractures and cracks 
in soil or porous rocks. One practical interest is that how long it would take for the pollutant 
to reach the whole (but finite) media or the typical saturation time of process. Obviously, the 
speed and efficiency of the diffusion process depends on the permeability, and the permeability 
in turn depends on the structure of porous rocks and the way of connections. The fraction or 
probability p of the cracks may have strong influence on the permeability even for the very small 
p <^ 1, which is a typical feature of small- world networks. 

The evolutionary solution of V{t) with time r and time delay 5 gives the volume affected by 
the pollutant. Solution (5) is valid for an infinite network. However, in reality most systems are 
finite. Since V{t) increases dramatically as r increases, we can expect that the whole system 
will be reached and saturated when t tends to a typical value Ts, which represents a system 
saturation time. For a finite system with a size N L"", the characteristic saturation time Ts 
can be estimated by 

Ts/S 

^ ^ {nj)\ ■ 

In most random walk models, it often the time delay is small or 5 — >■ 0. Now for simplicity, we 

take 5 = 0. For n = 1, we have 

r, ~ ln(L/e) (8) 
or in terms of real time t ~ {^/A) ln(L/^) where A = (Tdin — 1)!)^''". For n = 2, we have 

Ts^cosh-^L'/f), or ~4ein(L/0. (9) 

In general, r,, (or ts) is a function of typical length scale ^. Substituting the express ^ in terms 
of p into the above equation, we have 

t.^'-^^. (10) 

pl/n ^ ' 

This shows that the fraction p of cracks in disordered porous rocks can influence on the saturation 
time significantly. For a fixed large L, then ts ^ (lnL)/p^^". Thus, the increase of p, say, 
p = 10~^ to p = 10~^ will decrease ts by a factor of 10 for n = 1. This have important 
implications in management of the waste disposal and pollutant control. For example, the 
pollutant wastes were disposed in the area with possible cracks even the possibility is very 
small, the risk increases dramatically as the probability increases as seen from the saturation 
time. 

This model can also apply to other processes such as the forest fire, oil recovery and saltwater 
intrusion. If there is a fire start in a forest, the sparking effect will increase the probability 

of the long range shortcuts, and thus decrease the saturation time and the fire spread more 
rapidly. Similarly, saltwater intrusion proceed more quickly in the permeable zone with even 
small portion of cracks. 



Multifractals 

Since the velocity of a particle travels at a constant speed u = 1, we can rewrite the solution 
(5) in terms of the distance r and r = rA/^ with A = (Td,(n — 1)!)^^", we have 

Ar/(tS) 

Since V(r) increases as r increases or V{r) ~ . One can expect that fractal dimension of the 
small-world networks will give a measure of the network properties as first done by Newman 
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Figure 1: Fractal dimension of small- world networks with time-delay i5 = 0, 1,5, 10 for a network 
size N = 500,000, r = 1000 and n = 1. Numerical results (marked with o, o etc) agree well with 
analytical express (solid) calculated from (11). 

and Watts [1999] which does not include the effect of time-delay. Then the fractal dimension D 
can be calculated by differential formulae 



D ■- 



d\iiV{r) 
dlnr 



(12) 



which can be calculated using (11). We now show how the time-delay affect the fractal dimension 
D of the small-networks. For 5 Q and n = 1, we have Z) — >■ 1 for r/^ ^ 1 and D — >■ A£^~^ 
for r/^ — > oo. For 5 7^ 0, there is no general explicit asymptotic expressions and numerical 
evaluations can be easily done. The dependence of the fractal dimension D on the length scale 
^ and the time-delay 5 suggests the multifractal features of the small-networks. 

In order to compare the analytical solution with numerical simulations, we also use numerical 
simulations [ VFaite,1999]. The numerical simulations for n = 1 and the comparison of fractal 
dimension with the analytical solution (11) is shown in Figure 1 where the network size A'^ = 
500,000, r = 1000 and 5 = 0, 1,5, 10. We can see that numerical results (marked with o, o 
etc) are in good agreement with the analytical solution (solid lines) by using (11). As expected, 
for p ~ 1/C ^ 10~^ (or ^ r so that r/,f ^ 1), both numerical results and analytical fractal 
dimension approach to n = 1. The time-delay has a very strong effect on the fractal dimension 
D of the small- world networks. For 5 S> 1, -D is substantially reduced compared with the one 
without time-delay 5 = 0, and D is quite near n for most of the region. The large time-delay 
actually makes the network become larger, and the influence spread slowly in more localized 
areas. 



Percolation in Porous Media 

Diffusion and transport process can be considered as a percolation phenomenon. In addition, 
permeability and conductivity can also be modelled as a percolation process on the intercon- 
nected networks by using Monte Carlo simulations or fractional Brownian motion method. From 
the discussion in the previous sections, we can see naturally that the small-world network model 
can show newlight in the modelling of subsurface flow and contaminant transport in disordered 
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porous media [Holzbercher, 1998]. In small- world model, all the site or vertices are susceptible 
to pollutant particles. However, in reality, there is only a fraction of the sites can do so as in 
the case of site percolation. In other words, each site is occupied with a probability or density 
of (j)- Site percolation can also be used as the idealized model for the distribution of oil/gas 
or hydrocarbon inside porous rocks. The average concentration of fluids in the porous rocks is 
presented by the occupation probability i^, which happens to be the porosity in the saturated 
sediments or shales. A test particle or ant in the labyrinth move randomly on the intercon- 
nected network and various studies [e.g., Stauffer and Aharony, 1994; Yang, 2000] show that the 
percolation pattern has the fractal feature. Now we shall discuss how 4> affect the behaviour of 
the small-world model. 

With site occupation probability as an additional parameter, the solution (5) shall be 
modified as 

At/(f6) 

Vyt) = <p5 2^ > <P > 9c, (13) 

where (j)c is the critical probability and the above solution is only meaningful for ((> > (pc- The 
critical value <^c can be estimated numerically. However, for some type of lattice structure, it 
can be caculated exactly (e.g., for 1-D lattice, 0c = 1)- The introduction of a number oi pkN 
shortcuts (where A'^ is the network size) will generally change 0c as discovered by Newman and 
Watts [1999] by solving the 1-D case 

For 2-D or higher dimension case, there is no general explicit solution for (j)c and numerical 
approach is usually used. For a lattice in n-dimensional percolation structure, formation of one 
very large cluster gives 

0cPfeiV=^^fffc,<A*(l-0c)1, (15) 

s 

or ^ 

p = -^'^ 9ks(t>c~'^ {1- - <i>cy], (16) 

s 

where the summation is over s which is the cluster perimeter [Stauffer and Aharony, 1994]. The 
coefficient Qks depends on the network or lattice structure. For example, the triplet on a square 
lattice corresponds to = 3 and s = 8 with gts = 2 or s = 7 with gks = 4. A non-zero p usually 
gives a lower value of 4>c than that for p = 0, which means that cracks or higher permeability 
connection will reduce the critical site occupation probability and thus increase the efficiency of 
the transport and diffusion process in disordered porous media. 



Discussion 

A modified small- world network model has been presented here to characterise the effect of time- 
delay on small-world networks. Numerical simulations and analytical analysis for networks with 
time-delay show that the time-delay parameter 5 has a very strong effect on the fractal dimension 
and other properties such as the speed and saturation time of the delay networks. A small-world 
network can become larger if sufficient time-delay in the system response in introduced. 

On the other hand, in order to make a large world network transfer into a small-world or 
to increase the transport and diffusion efficiency of the network, a slight higher probability p of 
long-range random shortcuts are necessary compared with the one without time-delay because 
this essentially reduce the characteristic length scale ^ = 0{l/p) quite significantly. Cracks 
and other long-range connections can have a strong influence on permeability, conductivity and 
percolation in disordered porous media. Further studies of the dynamics of small- world networks 
are obviously necessary and the present work is just one of such attempts to the understanding 
of dynamics of networks and application in geophysical modelling. 
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